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Abstract: We examine the prospects for super symmetry discovery in the minimal super- 
gravity (mSUGRA) model via indirect detection of neutralino dark matter. We investigate 
rates for muon detection in neutrino telescopes, and detection of photons, positrons and 
anti-protons by balloon and space based experiments. We compare the discovery reach in 
these channels with the reach for direct detection of dark matter, and also with the reach 
of collider experiments such as Fermilab Tevatron, CERN LHC and a yfs = 0.5 — 1 TeV 
linear e + e~ collider. We pay particular attention to regions of model parameter space in 
accord with recent WMAP results on the dark matter density of the universe. We find 
that 3rd generation direct dark matter detection experiments should be able to cover the 
entire WMAP allowed portion of the hyperbolic branch/focus point (HB/FP) region of 
parameter space, while the IceCube neutrino telescope can cover almost all this region. 
This is in contrast to the case of the CERN LHC or a linear e + e~ collider, where only a 
fraction of the HB/FP region can be accessed. In addition, we show that detection of 7s, 
e + s and ps should occur in much of the HB/FP region, as well as in the low mi/ 2 portion 
of the A annihilation funnel, and will be complementary to searches via colliders in these 
regions. 
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1. Introduction 



Evidence for cold dark matter (CDM) in the universe comes from observations of galac- 
tic rotation curves and binding of galaxies in clusters, from matching observations of 
large scale structure with simulations, from gravitational microlensing, from the bary- 
onic density of the universe as determined by Big Bang nucleosynthesis, from observa- 
tions of supernovae in distant galaxies, and from measurements of anisotropics in the cos- 
mic microwave background radiation (CMB)(l|. In particular, the Wilkinson Microwave 
Anisotropy Probe (WMAP)|§] collaboration has extracted a variety of cosmological pa- 
rameters from fits to precision measurements of the CMB radiation. The properties of a 
flat universe in the standard AC DM cosmological model are characterized by the density 
of baryons (£lb), matter density (f2 m ), vacuum energy (^a) and the expansion rate (h) 
which are measured to be: 

n b = 0.044 ± 0.004 (1.1) 
n m = 0.27 ±0.04 (1.2) 
Q A = 0.73 ±0.04 (1.3) 

-0.04 
-0.03- 



h = 0.7ll°;°£. (1.4) 



From the WMAP results, a value for the cold dark matter density of the universe can be 
derived: 

VcDMh 2 = 0.1126lS(taM8i) at 68 ( 95 )% CL - (1-5) 
While the origin of dark energy in the universe remains a conundrum, there exists a number 
of hypothetical candidate elementary particles to fill the role of CDM. 

A particularly attractive candidate for CDM is the lightest neutralino in i?-parity 
conserving supersymmetric models Q. In the paradigm minimal supergravity (mSUGRA) 
model y], it is assumed that at the scale Q = Mqut-, there is a common scalar mass 
m,Q, a common gaugino mass m\ji, and a common trilinear term A§. The soft SUSY 
breaking terms can be calculated at scale Q = M wea }~ via renormalization group evolution. 
Electroweak symmetry breaking occurs radiatively (REWSB) due to the large top quark 
mass, so that the bilinear soft breaking term B can be traded for the weak scale ratio of 
Higgs vevs tan /3, and the magnitude (but not the sign) of the superpotential [i term can 
be specified. Thus, the mSUGRA model is characterized by four parameters plus a sign 
choice: 

mo, mi/2, A); tan/3, and sign(fi). (1-6) 

Once these model parameters are specified, then all sparticle masses and mixings are de- 
termined, and scattering cross sections may be reliably calculated. 

In the early universe at very high temperatures, the lightest neutralino Z\ will be in 
thermal equilibrium, so that its number density is well determined. As the universe expands 
and cools, the expansion rate outstrips the neutralino interaction rate, and a relic density 
of neutralinos is frozen out. The neutralino relic density ^^h 2 at the present time can be 
determined by solving the Boltzmann equation for neutralinos in a Friedmann-Robertson- 
Walker universe. 
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In most of the parameter space of the mSUGRA model, it turns out that a value of 
^Zi^ 2 we ^ beyond the WMAP bound is generated. Only certain regions of the mSUGRA 
model parameter space give rise to a relatively low value of ^^ ± h 2 in accord with cosmo- 
logical measurements and theory. These regions consist of: 

1. The bulk annihilation region at low values of tjiq and m^, where neutralino pair 
annihilation occurs at a large rate via i-channel slepton exchange. 

2. The stau co-annihilation region at low mo where — so that Z±s may co- 
annihilate with fis in the early universe ||. 

3. The hyperbolic branch/focus point (HB/FP) region|| at large mo near the boundary 
of the REWSB excluded region where |//| becomes small, and the neutralinos have 
a significant higgsino component, which facilitates annihilations to WW and ZZ 
pairs g |. 

4. The ^-annihilation funnel, which occurs at very large tan/? ~ 45 - 60 g. In this 
case, the value of m^ ~ 2m^. An exact equality of the mass relation isn't necessary, 
since the A width can be quite large (r^ ~ 10 — 50 GeV); then 2m^ can be several 
widths away from resonance, and still achieve a large Z\Z\ — > A — > ff annihilation 
cross section. The heavy scalar Higgs H also contributes to the annihilation cross 
section. 



In addition, there exists a region of neutralino top-squark co- annihilation |i0| (for very 
particular Aq values) and a light Higgs h annihilation funnelQ (at low m^ values). 

In past years, the bulk annihilation region of parameter space was favored. This sit- 
uation has changed in that the low mo and m-i/2 portion of the bulk annihilation region 
has been excluded by LEP2 chargino and Higgs search bounds, while the larger mo and 
m i/2 portion generally predicts values of Q^^ 2 beyond the rather restrictive upper bound 
of ^h 2 < 0.129 obtained from WMAP. Any remaining portions of the bulk region give 
rise to large- usually anomalous- predictions of the rate for BF(b — > S7) decays and muon 
anomalous magnetic moment = (g — 2)^/2[|ll], An increase of either of the pa- 
rameters mo or mi/2 leads generally to heavier sparticle masses and values, so that 
predictions for loop induced processes become more SM-like. 

A panoply of collider and non-accelerator experiments are now operating or will soon 
be deployed that will shed light on CDM. Prospects for detecting dark matter and deter- 
mining its properties are particularly bright in the case of neutralinos from supersymmetry. 



Neutralino dark matter may well be produced at the Fermilab Tevatron pp collider 1 13], the 



CERN LHCQ pp collider, and a yfs = 0.5 - 1 TeV linear e + e~ collider 0]. In addition, 
there exist both direct and indirect non-accelerator dark matter search experiments that 
are ongoing and proposed. Direct dark matter detection has been recently examined by 



many authors [16], and observable signal rates are generally found in either the bulk anni- 
hilation region, or in the HB/FP region, while direct detection of DM seems unlikely in 
the A- funnel or in the stau co-annihilation region. 

Indirect detection of neutralino dark matter [fT^| may occur via 
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1. observation of high energy neutrinos originating from Z\Z\ annihilations in the core 
of the sun or earth [18], 

2. observation of 7-rays originating from neutralino annihilation in the galactic core or 
halo HT9|1 and 



3. observation of positrons [^] or anti-protons [21] originating from neutralino annihila- 
tion in the galactic halo. 

The latter signals would typically be non-directional due to the influence of galactic mag- 
netic fields, unless the neutralino annihilations occur relatively close to earth in regions of 
clumpy dark matter. 

The indirect signals for SUSY dark matter have been investigated in a large number of 



papers, and computer codes which yield the various signal rates are available[22, 23]. Re- 
cent works find that the various indirect signals occur at large rates in the now disfavored 
bulk annihilation region, and also in the HB/FP region Q. Naively, this is not surprising 
since the same regions of parameter space that include large neutralino annihilation cross 
sections in the early universe should give large annihilation cross sections as sources of indi- 
rect signals for SUSY dark matter. In Ref. p4|] , it was pointed out that the A annihilation 
funnel can give rise to large rates for cosmic 7s, e + s and ps. However, neutralino-nucleon 
scattering cross sections are low in the A annihilation funnel, so that no signal is expected 
at neutrino telescopes, which depend more on the neutralino-nucleus scattering cross sec- 
tion than on the neutralino annihilation rates. Our goal in this paper is to combine the 
projected discovery contours from Tevatron, LHC and LC searches with those of direct 
and indirect dark matter search experiments. It turns out that each distinct region of 
mSUGRA parameter space which gives rise to an acceptable O^/t 2 value also gives a set 
of unique predictions for combinations of collider and non-accelerator experiments. 



We begin our analysis by generating sparticle mass spectra using Isajet v7.69|2q|, which 
includes full one-loop radiative corrections to all sparticle masses and Yukawa couplings, 
and minimizes the scalar potential using the renormalization group improved 1-loop effec- 
tive potential including all tadpole contributions, evaluated at an optimized scale choice 
which accounts for leading two loop terms. Good agreement between m/j values is found in 
comparison with the FeynHiggs program, and there is good agreement as well in the uia 
calculation between Isajet and SoftSUSY, Spheno and Suspect codes, as detailed in Ref. 



1 26]. To evaluate the indirect signals expected from the mSUGRA model, we adopt the 



DarkSUSY 3.14 package |23|] interfaced to Isasugra 1 . For our calculation of the neutralino 
relic density, we use the IsaReD program[27] interfaced with Isajet. IsaReD calculates all 
relevant neutralino pair annihilation and co-annihilation processes with relativistic ther- 
mal aver aging p^]. An important element of the calculation is that IsaReD calculates the 
neutralino relic density using the Isajet t, b and r Yukawa couplings evaluated at the scale 
Q = rriA- The Yukawa coupling calculation begins with the DB,' fermion masses at scale 
Q = Mz, and evolves via 1-loop SM renormalization group equations (RGEs) to the scale 



1 Isasugra is a subprogram of the Isajet package that calculates sparticle mass spectra and branching 
fractions for a variety of supersymmetric models 
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Qsusy = yJ m t L m t R j where complete MSSM 1-loop threshold corrections are implemented. 
Evolution at higher mass scales is implemented via 2-loop MSSM RGEs. The final RGE 
solution is gained after iterative running of couplings and soft terms between Mz and 
Mqut and back until a convergent solution is achieved. 

The remainder of this paper is organized as follows. In Sec. ||, we present an overview 
of indirect, direct and collider search experiments for neutralino dark matter. In Sec. ^, we 
discuss the impact of different halo models on our calculations. We then present our main 
results, which are a series of plots in the tuq vs. plane of the mSUGRA model showing 
the comparative reach of various indirect, direct and collider searches for neutralino dark 
matter. In Sec. [|, we present our conclusions. 



2. Experimental overview 
2.1 Neutrino telescopes 

A novel technique for detecting dark matter is to search for neutrino signals coming from 
neutralino annihilation in the core of the earth or the sun @. As the sun or earth proceeds 



on its orbital path, neutralinos can be swept up and become gravitationally captured by 
the process of neutralino-nucleon scattering until the recoil neutralino velocity drops below 
the escape velocity. Thus, a high density of neutralinos can accumulate in the core of the 
earth or sun, where they can efficiently annihilate. Neutralino annihilation into SM final 
states such as bb, cc, W + W~~ or ZZ ultimately yields neutrinos via the b — > c£u£, c — ► s£ve, 
W — > Iv or Z — > vipn decays. The neutrinos can propagate out of the core of the earth 
or sun, and be detected via — ► \i conversions in neutrino telescopes such as Antares or 
IceCube. In fact, limits have already been obtained by Amanda for the case of neutralino 
annihilation in the core of the earth [p9|]. 

The Antares v telescope should be sensitive to En > 10 GeV; it is in the process 



of deployment and is expected to turn on in 2006 [3C]. It should attain a sensitivity of 
100 — 1000 /is/km 2 /yr. The IceCube v telescope is also in the process of deployment at 
the south pole^TJ. It should be sensitive to Ep > 25 — 50 GeV, and is expected to attain 
a sensitivity of 40 — 50 /is/km 2 /yr. Full deployment of all detector elements is expected to 
be completed by 2010. 

The rate for neutralino annihilation in the sun or earth is given by 

r A = ^tanh 2 (t /r), (2.1) 

where C is the capture rate, A is the total annihilation rate times relative velocity per 
volume, t® is the present age of the solar system and r = 1/V CA is the equilibration time. 
For the sun, the age of the solar system exceeds the equilibration time, so Ta ~ and the 
muon flux tends to follow the neutralino-nucleon scattering rate rather than the neutralino 
pair annihilation cross section. In this case, the indirect dark matter detection rate should 
be relatively independent of the dark matter halo profile, aside from the value of the local 
neutralino relic density. Thus, uncertainties in the predicted rates for neutrino detection 
via neutralino annihilation in the core of the sun should be low. In contrast to the sun, the 
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earth has typically a much longer equilibration time r, so that Ta ~ \C 2 At 2 , and is hence 
more sensitive to the neutralino annihilation cross section times relative velocity. Muon 
rates from neutralino annihilation in the core of the earth are typically much lower than 
those from the sun. In addition, expected rates from the earth may be diminished even 
further by solar depletion effects: see J. Lundberg and J. Edsjd, Ref. |3^] (these effects are 
not included in DarkSUSY 3.14, but are included in DarkSUSY 4.0). 

2.2 Detection of 7s 

Neutralinos may also collect in the core of the galaxy, where they can annihilate at a high 
rate. In this case, Z\Z\ — > qq, W + W~, ZZ — ► hadrons which gives rise to photons 
typically from ir° — > 77 decay. It is also possible for Z\Z\ — ► 77 (or Z-y), in which case 
E~f ~ m Z\ ^ e s ig na ture is spectacular in this case, but the rates are loop suppressed. 
The 7 rays can be detected down to sub-GeV energies with space-based detectors such 



as EGRET [33 1 or GLASTJ3J]. Ground based arrays require much higher photon energy 



thresholds of order 20 — 100 GeV. Experiments such as GLAST should be sensitive to rates 
of order 10 -10 7s/cm 2 /sec assuming E~ > 1 GeV. In fact, it has recently been suggested 
that the extra-galactic gamma ray background radiation as measured by EGRET is well 
fit by a model of neutralino annihilation^]; see also |3(|. It is important to note that the 
prediction for rates for 7 detection depends sensitively on models for the neutralino density 
near the galactic core. The latter quantity is poorly known, so there can be a wide range 
in predicted rates depending on assumptions about the galactic halo profile. 

2.3 Detection of e + s 

Cosmic positrons may also be searched for from neutralino annihilations in the galactic 
halo. In this case, the positrons would arise as decay products of heavy quarks, leptons and 
gauge bosons produced in neutralino annihilations. Space based anti-matter detectors such 



as Pamela[ 37] and AMS-02[38] will be able to search for anomalous positron production 



from dark matter annihilation. The cosmic positron excess as measured by HEAT [39] 



has been suggested as having a source in galactic halo neutralino annihilations [40, 36]. 
It is suggested by Feng et al. 0] that a reasonable observability criteria is that signal-to- 
background (S/B) rates should be greater than the 1 — 2% level. To calculate the S/B 
rates, we adopt fit C from Ref. (7j for the E 2 dQ e + /dCldE background rate: 

E 2 d$ e +/dndE = 1.6 x 10~ 3 E~ 123 , (2.2) 

where E is in GeV. We compute the signal using the DarkSUSY positron flux evaluated 
at an "optimized" energy of E = /2, as suggested in Ref. 0]. A S/B ~ 0.01 rate may 
be detectable |7|, 17 1 by experiments such as Pamela and AMS-02. 

2.4 Detection of ps 

Antiprotons may also be produced in the debris of neutralino annihilations in the galactic 
halo. Such antiprotons have been measured by the BESS collaboration^]. The differential 
flux of antiprotons from the galactic halo, d&p/ 'dEpdVL, as measured by BESS, has a peak in 
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the kinetic energy distribution at Ep ~ 1.76 GeV. The height of the peak at Ep ~ 1.76 GeV 
is ~ 2 x 10 -6 p/GeV/cm 2 /s/sr. Signal rates in the range of 1CP 7 — 1CP 6 p/GeV/cm 2 /s/sr 
might thus provide a benchmark for observability. 

2.5 Direct search for neutralino DM 

If indeed all space is filled with relic neutralinos, then it may be possible to directly detect 
them via their scattering from nuclei. Early limits on the spin-independent neutralino- 
nucleon cross-section (a SI ) have been obtained by the CDMS01, EDELWEISS @ and 



ZEPLIN1[44] groups, while a signal was claimed by the DAMA collaboration [f45|] . Collec- 
tively, we will refer to the reach from these groups as the "Stage 1" dark matter search. 
Depending on the neutralino mass, the combined limit on asi varies from 1CT 5 to 1CT 6 pb. 
This cross section range is beyond the predicted levels from most super symmetric mod- 



els. However, experiments in the near future like CDMS2, CRESST2p6f, ZEPLIN2 and 
EDELWEISS2 (Stage 2 detectors) should have a reach of the order of 10~ 8 pb. In fact, the 
first results from CDMS2 have recently appeared, and yield a considerable improvement 



over the above mentioned Stage 1 results [47]. Finally, a number of experiments such as 



GENIUS HI, ZEPLIN4j49| and XENON g] are in the planning stage. We refer to these as 
Stage 3 detectors, which promise impressive limits of the order of asi < 10 -9 - 10 -10 pb, 
and would allow the exploration of a considerable part of parameter space of many su- 
persymmetric models. In particular, the Stage 3 direct DM detectors should be able to 
probe almost the entire HB/FP region of mSUGRA model parameter space. We note 
here in addition that the Warm Argon Program (WARP)[[5l|l promotes a goal of detecting 
neutralino-nucleus scattering cross sections as low 1CP 11 pb. 

2.6 Fermilab Tevatron 

The Fermilab Tevatron pp collider can search for neutralino dark matter by detecting 
superparticle production reactions which lead to anomalous missing energy in the final 
state. For mSUGRA model parameter choices in accord with bounds on the chargino mass 
from LEP2 (fTt^ > 103.5 GeV), the most promising discovery channel is the clean trilepton 



plus Et final state [52], which typically originates from pp — * W\Z^X — » 3i+ Et + X, 
where X stands for assorted hadronic debris, and where the trileptons tend to originate 
from W\ — ► tvi>Z\ and Z2 — > UZ\ decay. Trilepton signal rates and backgrounds have been 
calculated in Ref. 53] using cuts SC2. These results were updated and extended to large 
mo regions of parameter space in Ref. [54]. 

2.7 CERN LHC 

The CERN LHC can search for neutralino dark matter by detecting superparticle pro- 
duction reactions which lead to anomalous missing energy in the final state. For LHC, 
however, gluino and squark production reactions are expected to be the dominant SUSY 
cross sections. The gluinos and squarks can decay through possibly lengthy cascade decays 
so that signal events will consist of multi-jets plus isolated leptons plus -^r|55|]. Signal and 
background levels have been computed in Ref. [pfl] for various combinations of jet and 



-6- 



lepton cuts. We adopt the most recent calculations of Baer et al. ||5q| for our projections of 
the LHC reach, assuming 100 fb _1 of integrated luminosity. 

2.8 Linear e + e collider 

A linear e + e~ collider operating at CM energy y/s = 0.5 — 1 TeV can also search for neu- 
tralino dark matter by detecting superparticle production reactions which lead to anoma- 
lous missing energy in the final state. The ultimate reach limits in the mSUGRA model 
depend on which sparticles are kinematically accessible to production. The reach contours 
are determined by[57, [5^] 



• e + e — * £ + £ at low mo, followed typically by I — » IZ\ decay and 

• e + e~ — > Wi~Wi production at moderate to high mo values, where typically W\ — > 
ff'Zi or WZi (/ is any SM fermion). 

• In intermediate regions, e + e~ — * Z1Z2 and/or e + e~ — * ZH, Ah may be accessible. 

In the HB/FP region, the superpotential n parameter becomes small, and the W\ 
and Z\ become increasingly higgsino like, and nearly mass degenerate. In the small mass 
gap region, conventional cuts oriented towards a substantial ~ m z 1 mass § a P mus t 
be replaced by new cuts. Ultimately, an e + e~ LC should be able to see chargino pairs 
essentially up to the kinematic limit for their production, over all of mSUGRA parameter 
space. Reach plots have been presented in Ref. |5q] assuming 100 fb _1 of integrated 
luminosity for both the y/s = 0.5 TeV and 1 TeV machines. 

3. Direct, indirect and collider searches in the mSUGRA model 

In this section, we evaluate the reach of various indirect, direct and collider searches for 
neutralino dark matter in the mSUGRA model. For our predictions of indirect neutralino 



detection rates, we use the DarkSUSY program|22], modified to interface with Isajet v7.69. 
The Isajet subprogram Isasugra is used to predict the sparticle mass spectrum and decay 
widths for various super symmetric models, including mSUGRA. 

3.1 Calculational overview and dependence on halo model 

While predictions for the reach of colliders for neutralino dark matter are firmly grounded in 
perturbative quantum field theory, many of the predictions for direct and indirect detection 
of dark matter depend on the assumed density profile of the galactic halo. 

For halo model dependence in the distribution of dark matter, we adopt the default 
DarkSUSY value: a spherically symmetric isothermal distribution given by 

P(r) = Po (r/r ° r V 7 (l + (ro/aD^ (3.1) 
(l + (r/a) a ) — 

where (a, f3, 7) = (2, 2, 0), ro = 8.5 kpc is the distance of earth to the galactic center, 
po = 0.3 GeV/cm 3 is the local dark matter density and a = 3.5 kpc is a distance scale. A 
number of other halo profiles are available in DarkSUSY, including 
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Navarro, Prenk, White profile with a = 20 kpc, with (a, /3, 7) = (1, 3, 1) |59(| , 
Moore et al. profile with a = 28 kpc, with (a, j3, 7) = (1.5, 3, 1.5) |Kl and 



• Kravtsov et al. profile with po = 0.6 GeV/cm 3 , a = 10 kpc, with (a, /3, 7) = 
(2, 3, 0.4) ||. 

The alternative distributions all use tq = 8.0 kpc as the sun galactocentric distance value. 

In Fig. g, we show in frame a) the various halo model distributions. While the distribu- 
tions are in close accord at r = 8.5 kpc (the distance of the sun to the galactic center), the 
profiles disagree strongly as r —* 0, reflecting our relative ignorance of the density of dark 
matter expected near the galactic center. In frame b), we show the neutralino relic density 
versus mSUGRA parameter mo, with mi/2 = 550 GeV, Aq = 0, fi < and tan/? = 50. The 
horizontal line shows the WMAP upper bound on In frame c), we show the flux of 

muons with E^ > 25 GeV originating from neutralino annihilation to neutrinos in the core 
of the sun. These curves only depend on the local neutralino relic density, and not on the 
global galactic halo profile. Thus, the first three halo profiles yield the same predictions for 
the muon detection rate, while the Kravtsov distribution, with a higher value for the local 
relic density, also gives higher rates for neutrino telescopes. In frame d), we show the flux 
of photons <J> 7 (7s/cm 2 /sec) with E~ > 1 GeV emanating from the galactic center, within 
a cone of solid angle 0.001 sr. The rates are large at large mo in the HB/FP regionQ, 
and are also large at mo ~ 1000 GeV, where neutralino annihilation in the galactic core 
can occur through the A and H resonances [p^j. Overall, the predicted rates at fixed mo 
vary over several orders of magnitude, reflecting the different model predictions for the 
DM density at the galactic center. In frame e) , we show as well the expected positron 
signal-to-background (S/B) rate for positrons originating in neutralino annihilations in the 
galactic halo. We see that observable rates may again occur in the HB/FP region, and also 
in the ^4-annihilation funnel. However, the uncertainty in the predictions due to variation 
in the halo model is less severe than in the 7 case, since now the positrons are expected to 
arise relatively nearby in the galaxy, where the DM density is much more constrained. In 
frame j), we show the differential flux of antiprotons from the galactic halo, d<&p/dEpd$l, 
for Ep = 1.76 GeV, which is near the location of the peak in the p kinetic energy distribu- 



tion as measured by the BESS collaboration! 41]. We see from the figure that the largest 



rates occur in the HB/FP region, and also in the A-annihilation funnel; in these regions, 
the signal rates can extend into the region of observability. Again, the predicted rates 
are sensitive to the assumed halo model. We note, however, that the DarkSUSY default 
halo distribution tends to give the most conservative of rate predictions, and it is possible 
that indirect detection rates could be much higher than those shown using the default halo 
model. 

One final qualifying note: it is common practice to rescale direct or indirect detection 
rates by a factor Q^^ti 2 / (QcDMh 2 ) re f when ^^/i 2 < (&>CDMh 2 ) re f, where (QcDMh 2 ) re f is 
some reference value ~ 0.025 — 0.1 which would give rise to the assumed local dark matter 
density p = 0.3 GeV/cm 3 . We do not apply this practice here. It is possible in a variety 
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of non-standard cosmological models (such as those containing quintessential scalar fields 
or primordial anisotropics) to obtain a dark matter relic density much higher than that 
obtained in a KCDM model. See e.g. Ref. |3|] for details. 

3.2 Results for the mSUGRA model 

Our first results for the mSUGRA model are shown in Fig. ^ in the too vs. m^j plane 
for tan/3 = 10, Aq = and /i > 0. We take rat = 175 GeV. The left-most red region is 
excluded because the stau becomes the LSP (in violation of search limits for stable charged 
or colored relics from the Big Bang), while the right-most red region is excluded due to a 
lack of REWSB. The lower yellow region is excluded by LEP2 searches for chargino pair 
production (m^ > 103.5 GeV), while the region below the yellow contour is excluded by 
LEP2 Higgs searches (ra^ > 114.4 GeV for a SM-like Higgs boson). The green shaded 
region has Q^h? < 0.129, in accord with the upper bound on CDM from WMAP. The 
left-most green strip along the low mo excluded region is the stau co- annihilation corridor, 
while the right-most green region corresponds to the HB/FP region. A remaining green 
region sitting just atop the LEP2 excluded region is the light Higgs annihilation corridor, 



where 2m^ ~ 77^(631. The uncolored regions all give too large a CDM relic density, and 
are thus excluded. 

The Fermilab Tevatron reach contour corresponding to a 5<r sig nal for 10 fb" 1 of 
integrated luminosity is denoted by TEV, while the CERN LHC 5<r reach for 100 fb _1 
is denoted LHC. The reach of a yfs = 500 GeV and 1000 GeV linear collider is denoted 
by LC500 and LC1000, respectively, assuming 100 ftr 1 for each. We see from the figure 
that the Tevatron can cover the light Higgs annihilation corridor, while the LC1000 and 
LHC can cover the stau co-annihilation region. The LHC can cover the HB/FP region 
up to rrii/2 ~ 700 GeV, which corresponds to a reach in ra g of about 1.8 TeV. In this 
region, squarks and sleptons are 4-7 TeV in mass, and effectively decoupled from LHC 
production. LHC can generate an observable signal cross section provided m g is light 
enough. If ra g is too heavy (~ 1.8 TeV), then pp — > ggX occurs at too low of a rate. 
Charginos and neutralinos can be quite light in the HB/FP region since |//| is small, but 
their signal events are difficult for LHC to extract from background, owing in part to a 
decreasing ~ m z 1 mass § a P as IH decreases. The LC500 and LC1000 can still generate 
Wi Wi pairs in the HB/FP region if they are kinematically accessible, and should be able 



to extract the corresponding low energy release events above SM background 1 58]. In the 
HB/FP region, we have the unusual situation that the LC reach exceeds that of the CERN 
LHC. However, the high m.]y 2 portion of the HB/FP region gives rise to cases where the 
chargino mass is too heavy to be accessible by a ^fs = 1 TeV LC, so a thorough search for 
SUSY by colliders in this DM allowed region apparently can't be made (unless a higher 
energy LC is constructed). 
We also show contours of 

• Stage 3 direct detection experiments (a si > 10~ 9 pb; black contour), 

• reach of IceCube v telescope with & sun (/i) = 40 /is/km 3 /yr and > 25 GeV 
(magenta contour), 
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• the $(7) = 10 10 7s/cm 2 /s contour with EL > 1 GeV in a cone of 0.001 sr directed 
at the galactic center (dark blue contour), 

• the S/B > 0.01 contour for halo produced positrons (blue-green contour) and 

• the antiproton flux rate 3>(p) = 3 x 10~ 7 j?s/cm 2 /s/sr (lavender contour). 

As noted by Feng et aZ.Q, all these indirect signals are visible inside some portion of the 
HB/FP region, while none are visible in generic DM disallowed regions (under the assumed 
smooth halo profiles). The intriguing point is that almost the entire HB/FP region (up 
to mi/2 ~ 1400 GeV) can be explored by the cubic km scale IceCube v telescope! It can 
also be explored (apparently at later times) by the Stage 3 direct DM detectors. Given 
the relative time scales of the various search experiments, if SUSY lies within the upper 
HB/FP region, then it may well be discovered first by IceCube (and possibly Antares), 
with a signal being later confirmed by direct DM detection and possibly the TeV scale 
linear e + e~ collider. There is also some chance to obtain indirect 7, e + and p signals 
in this region. Notice that if instead SUSY lies within the stau co-annihilation corridor, 
then it will be easily discovered by the LHC (for tan/3 = 10), but all indirect detection 
experiments will find null results in their DM searches. 

In Fig. |3|, we show again the tuq vs. mi/2 plane, but this time for tan/3 = 30 (other 
parameters remain the same). The results are rather similar to the tan/3 = 10 case shown 
in Fig. ||, although the stau annihilation corridor has expanded somewhat to higher mj/j 
values, and also the low mo and (bulk) region of parameter space has become more 

accessible to direct DM searches and even indirect DM searches in the p channel. The 
CERN LHC can still cover the entire stau co-annihilation region, and the low mi/2 portion 
of the HB/FP region. The LC1000 can cover much of the stau co-annihilation region and 
much of the high HB/FP region. However, a search of the m^ft < 1400 GeV portion 
of the HB/FP region can be made by the IceCube neutrino telescope, and later, Stage 3 
direct DM search experiments can cover the entire region. 

Fig. |] shows the mo vs. plane for an even higher tan/3 = 52 value. In this 

case, the stau co-annihilation corridor has increased to well beyond the LHC reach, and in 
fact this region of mSUGRA parameter space appears to be one which is consistent with 
WMAP relic density bounds, but beyond reach of any of the planned collider, direct or 
indirect search experiments for neutralino dark matter. At this high of a tan/3 value, the 
value of pseudoscalar Higgs mass tha has dropped |34j to such a level that there exists a 
large amplitude for off mass shell Z\Z\ — > A* — > bb annihilation in the early universe 0]. 
This additional amplitude helps to amplify the low mo DM allowed region, where now a 
combination of slepton co-annihilation, off shell A, H resonance annihilation, and i-channel 
neutralino annihilation via relatively light sfermions all serve to expand the DM allowed 
region. However, the HB/FP region is relatively insensitive to changes in tan/3, and is 
still well covered by ^-telescopes and direct DM searches. The reach of direct DM search 
experiments has vastly increased in the low mo and my 2 regions, where the reach of Stage 
3 detectors is comparable to the LC1000. The increase in direct DM detection rates is due 
to enhanced scattering via t-channel Higgs exchange graphs; these are proportional to the 
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square of b and r Yukawa couplings, which are large at large tan/?||. In addition, at low 
m,Q and mi/2, the DM allowed region may give rise to detectable rates for 7, p and e + 
detection. 

Fig. [| shows the uiq vs. mu2 plane for tan/3 = 55. In this case, the A annihilation 
funnel has moved into the central part of parameter space, opening up a large new region 
that gives a neutralino relic density in accord with WMAP results ||. The A annihilation 
funnel extends beyond the 100 fb _1 reach of the LHC, and makes a case for LHC running 
with much higher integrated luminosities, if no physics beyond the SM is found. In addition, 
the searches for 7s, e + s and ps are all enhanced in ths region, since now halo neutralinos 



can also annihilate through the broad A and H resonances 1 24]. Searches in these channels, 
however, cover only the low m^j portion of the A annihilation funnel. Annihilation through 
the A funnel also serves to expand somewhat the breadth of the HB/FP region. Even so, 
most of the HB/FP region can still be covered by the IceCube v telescope, while Stage 3 
DM detectors can cover the entire region. Much of it can also be covered by LHC, LCs and 
the search for halo annihilations into 7s, e + s and ps. As tan j3 increases much beyond 55, 
the parameter space starts to collapse due to inappropriate breakdown of EW symmetry. 

In Fig. |(| we show the mo vs. ra^ plane for fi < 0, Aq = and tan/3 = 45. For 
negative \i values, the ^-annihilation funnel enters the plane at lower tan/3 values; it is 
hence narrower than in Fig. [|, since the b and r Yukawa couplings are smaller, and the 
A and H widths not so wide. In this case, the CERN LHC covers almost all the stau co- 
annihilation region and the A funnel, and would certainly cover all of these with a higher 
integrated luminosity. As before, the HB/FP region is only partially covered by LHC and 
also by a LC, but again it is covered by the IceCube v telescope up to mx/2 ~ 1400 GeV, 
and covered completely by Stage 3 direct DM detectors. In addition, the enhanced rates 
for 7s, e + s and ps are displayed in the A annihilation funnel. 

We show in Fig. ^ the same ijiq vs. plane for {i < 0, but this time for tan j3 = 50. 
In this case, the A annihilation funnel is more centrally located, and shows observable 
rates for 7s, e + s and ps for the lower portion of the funnel. The red bulge moving into the 
figure at low mo and low mu2 denotes the region where m\ < 0, and begins the collapse 
of parameter space at high tan/3. The HB/FP region is again almost completely covered 
by IceCube, although the low m 1 / 2 portion of this region is no longer accessible to direct 
DM searches, owing to interferences in the neutralino-proton scattering rates. 

4. Conclusions 

In previous reports, the reach of the Fermilab Tevatron, the CERN LHC and y/s = 0.5 
and 1 TeV e + e~ linear colliders has been computed in the mSUGRA model. In addition, 
the reach of direct DM detection experiments has also been worked out, and found to be 



in many respects complementary to collider searches [ 35 1. In this paper, we augment these 
previous works by presenting as well the reach of various indirect search experiments for 
neutralino dark matter. These include searches for neutralino annihilation in the core of 
the sun (or earth), leading to detection of — ► /i conversions in neutrino telescopes such 
as Antares and IceCube. Also, we show reach contours for indirect searches for neutralino 



- 11 - 



annihilation in the galactic core via 7 detection, and for neutralino annihilation in the 
galactic halo via e + and p detection. Our results have focussed mainly on the WMAP 
allowed regions of the paradigm mSUGRA model. 
We have several main conclusions: 

1. In the stau co- annihilation region, indirect searches for neutralino dark matter have 
only a feeble reach. The best reach comes if tan /3 is large and there is some overlap 
with the A annihilation funnel. However, the CERN LHC can probe all the stau 
coannihilation region for tan /3 ~ 45. Much of it can also be explored by linear e + e 
colliders. In fact, the large mj/2 portion of this region at large tan/3 seems to be 
a region of mSUGRA parameter space which is consistent with WMAP limits on 
Jl^/i 2 , but not accessible to any planned experiments. 

2. Most of the A annihilation funnel can be explored by the CERN LHC, although again 
the large portion of it might not be accessible to any search experiments. The 
lower portion of the A funnel is also accessible to 7, e + and p searches for neutralino 
annihilation in the galactic core or halo. The indirect detection reach in this region 
is comparable to that of a y/s = 1 TeV linear e + e~ collider. Detection of v M — > \i in 
neutrino telescopes is unlikely to occur in this region. 

3. In the HB/FP region, the CERN LHC can cover m^ 2 values up to ~ 700 GeV, 
corresponding to a value of nig ~ 1.8 TeV. Linear e + e~ collider can do better, since 
they can detect chargino pair production, even if the energy release from chargino 
3-body decay is very low. The LC reach is limited by their CM energy, and a 1 TeV 
linear collider will not be quite sufficient to explore the entire high mi/2 portion of 
the HB/FP region. However, in this region, rates for detection of neutrinos arising 
from neutralino annihilation in the core of the sun are large, and it seems likely 
that IceCube can explore or rule out the HB/FP region for m\/2 < 1400 GeV. In 
addition, Stage 3 direct DM detection experiments sensitive to neutralino proton 
spin-independent scattering cross sections of 10 -9 pb should be able to access the 
entire HB/FP region, unless fi < and tan/3 is large, in which case a small hole 
arises in the low portion (which will be explored by LHC and LC anyway). 

Ultimately, the search for neutralino dark matter can proceed via direct DM searches, 
indirect DM searches and collider searches. By combining results, all the different search 
experiments can cover almost all the WMAP allowed mSUGRA model parameter space, 
save for a few regions which occur in the high mi/2 portion of the stau coannihilation 
corridor or the A annihilation funnel. 

Note added: A similar study of indirect and direct signals for neutralino dark matter 
in the mSUGRA model appeared shortly before the release of this work by Edsjo, Schelke 
and Ullio |]66|| . Where the two papers overlap, we seem to be in agreement. 
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Figure 1: In frame a), we show various halo model density profiles versus galactic radial distance 
r (kpc), while in frame 6), we show the neutralino relic density, along with the WMAP bound 
(horizontal line). Frame c) shows the flux of muons from the sun. In frames d), e), and f) we show 
rates for detection of 7s, e + s and ps from neutralino annihilations in the galactic core and halo, 
for different halo model choices. Frames b)-f) are plotted versus mo for mwa = 550 GeV, Aq = 0, 
tan j3 = 50 and [i < 0. 
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Figure 2: A plot of the reach of direct, indirect and collider searches for neutralino dark matter 
in the nig vs. m 1 / 2 plane, for Aq = 0, tan /3 = 10 and fi > 0. 
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Figure 3: A plot of the reach of direct, indirect and collider searches for neutralino dark matter 
in the toq vs. m 1 / 2 plane, for Aq = 0, tan [3 = 30 and fi > 0. 
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Figure 4: A plot of the reach of direct, indirect and collider searches for neutralino dark matter 
in the nig vs. m 1 / 2 plane, for Aq = 0, tan [3 = 52 and fi > 0. 
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Figure 5: A plot of the reach of direct, indirect and collider searches for neutralino dark matter 
in the nig vs. m 1 / 2 plane, for Aq = 0, tan [3 — 55 and fi > 0. 
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Figure 6: A plot of the reach of direct, indirect and collider searches for neutralino dark matter 
in the toq vs. m 1 / 2 plane, for Aq = 0, tan [3 — 45 and \l < 0. 
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Figure 7: A plot of the reach of direct, indirect and collider searches for neutralino dark matter 
in the mo vs. m 1 / 2 plane, for Aq = 0, tan [3 = 50 and fi < 0. 
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